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On the Values of the Constants in the Equation r A r x (r) + r Ar-i xCr 1 
'+ . . . + r A^ (0 + .. . 4 >A 0 —y x = o ; obtained by the method 
of least squares , from the n+ i values of y x ivhen x = o, 1, 2, 
, . . n ; n being greater than r. By C. Carpmael, Esq. 

Suppose it known, or assumed, that the value of a quantity 
y x (the approximate value of which is known when x = o, and 
when x — 1, 2, &c. . . . n) may be expressed as the sum of a 
series of factorials of the form r A t x lJ \ where ,.A t is the coefficient 
of 03® when the highest factorial in the series is 03 (r) . 

The object of the following investigation is to determine the 
values of the coefficients, r A r &c., -which make the sum of the 
squares of the differences between the given values of y and 
those obtained from the formula a minimum.* 

Now in order that 


A=0 + rK-l*< r ~» + • • • + + • • • + rh Q -y x f 


may be a minimum, r A., &c. must satisfy simultaneously the r -f 1 
equations of the form 


+ r A,_i S^r-n^ + . . . + r A,^U t . • • 


+ r hoKJo* (s) -KJo* (s) y* = o, 


s having in succession the values r, r— 1, r— 2, . . . , o. 

The value of r A r obtained from these equations is a fraction, 
whose denominator is a determinant of the r + ith order, having 


ar=n 


^(s) X (t) 


* The notation is : x^ = x ( x~ 1) . . . (x — r+ 1); \r = 
P* - Vn + Vn- 1 • . • +2/0.— Ed. 


.2.3. 
p p 


r ; 
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!for the constituent in its r — s-fi 1 th row and r—t -\-1 th column ; 
land whose numerator may be obtained from the denominator 
!by substituting constituents of the form 


- for those of the form 


x=n 

0 


x& y x 


f- n «x& xW 


in the first column. 

Performing the summations wherever y is not involved, by aid 
of the formula 


A" 1 u x v x = u c A- 1 v x - A u x A“ 2 v x+l + A' 2 u x A" 3 v x+2 - &e'., 


we obtain for the value of the constituent in the r — s + 1 th row 
and r—t+ 1 th column 


» +1 


it) n + 1 


0 +D 


s + 1 


■—t n + 1 


... • ##+ 2O+2) ] t \ s 

- + . . . + (-!)* - 3 ^— 


S+l s + 2 


+ i + 1 


-j- / (AS + t + 1), 



In the last column t = o, and (i) reduces to 


n + 1 


(s+l) 


8 + I 


in the last but one i = 1, and (i) becomes 


_ n+ i (s+I) n + z^ 

n + 1 -— -.• - : 

s+1 U1 8+2 


and so on. Subtracting ?i+i times the last column from the 
last but one, we obtain new determinants of the same values as 
the previous ones, and having the constituents in the last column 
but one of the form 


S + I H2 


Similarly, by subtracting 2 , n 4 -1 times this new column, and 
n -f 1 ,n times the last column from the last column but two, we 
reduce its constituents to the form 


2 


n + 


“Cx+3) 

.■) 


i -f I $- + 2 + 3 


By proceeding in this way we reduce the determinants, without 
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'CTl 1 

loo 1 

1^1 .. 

^hanging tbeir valu es, to ot hers having the constituents in the 
!^!~s+ 1 th row and r— t -f 1 th column of the form 



■[. 1 

loo 1 

I'—11 


(— - - - j - —n + t + i (l,+m) , 
|.s + t+ I - 


s and t having any value from o to r, except in the first column 
of the numerator, which is left as before. 

ISTow divide the r— s+ 1 th row in each of these new determi¬ 
nants by 


and the r — t + 1 th column by 

(_ l)t '1 n+t+i rt 

\r + t+ 1 


This wil l reduce the constituent in the r—s 4-1 th row and the 
r—t -f 1 th column to the form 

r + t+ I _ . 

7-,— " or !’ + rn v A "', 

} s + t + 1 ’ 

except in the first column of the numerator, where the r —5 + 1 th 
constituent will be 

\ 2 r+1 ^lZ n y x 

f- r )r h= - r=a J 

I L Li I+F*T°’ +m) 


As we here divide both the determinants by the same quantities, 
the value of their ratio is unchanged. 

The value of r A r thus obtained may be expanded in the form 


-A r = 2,“ 0 (-1) 


.. ! 2 r - i 


f *’ n + r - l : " 




(h) 


where r R a is of the form 


r p 2 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at University Library on March 16, 2015 



























xxxix. 9 ? 


392 

'CSV 

loo 1 
1^1 

■cTl 1 

lOO I 


' 4 *^ 

«c 2 ' 
1S1 

.2. 


M + + 

5 - S> 


Mr. Garpnael , On the Values of the 


y 




+ + 
5 . 




the numerator being the same as the denominator, except that 
the first column and r-s-t i th row are wanting. 

Now reverse the order of the columns, to do which will 
require \r r— 1 interchanges of adjacent columns in the nume¬ 
rator and -|r r-f-i in the denominator, and will therefore alter 
the sign of the fraction, or not, according as ( —i) r is positive or 
negative. 
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loo 1 
1^1 


Hence 


1S1 

,Si 


M S- 
N 


5 s 

W 


+ 

? 


Ss 


+ 

5s 


+ 


+ + 

5 . 5 . 


+ 

ot 


+ 

5s 


+ 

5 > 


M + 

5 » 


+ 

ft 


+ 

ft 


+ 

5 s 


I 


In the determinants in (iii) from each column, beginning at 
the right, subtract that next to its left; the determinants will 
thus be each reduced one order lower, for the constituents of the 
first row are all reduced to zero, except the first, ’which remains 
unity. Then divide each row by the factor common to that row. 
These common factors are the same in both numerator and de¬ 
nominator, except that one factor r — s occurs in the denominator 
and not in the numerator. The value of r R s is now reduced to 
the form 
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( 1 

r 4-1 


. 2T — I 


(r + 1 ) 0 —»-s) . . . 

( r+ l)0--^) . . . . 


. . (2i*-l)( r -*- 2 ) 

. . (2?’— I)(»—') 


(-0 r 

r — s 


( r+ l)(r-l) 


. . . . ( 2 r —^ 


/ I 

r + I 


I 

2 /’ 


(r + 1 i) 


(2r) 


\ (r + i)( r -i) . . , ..^ 

v^here the determinants are similar in form to those in (iii) but 
one order lower, and the missing row in the numerator is now 
that containing factorials of r— 5—1 factors. 

Repeat the reduction in this way r — s times, including that 
already performed, then the value of r R s will be reduced to the form 


' r +1 r + 2 ....... r + s 

(r + i)(*> (r + 2)®.(r + d)® 


. (r+i)® (r+2)® . 


r — s 

- - ( I I 


r+i r + 2 . 

(r + 1)® (r + 2)® 


(r + 5 )® ) 


r + s + I 
(r + s + 1 )® 


\ (r+l)® (r + 2)® ...... (r + s+l)® 
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■(Ti 1 

loo 1 

1^1 _ 

! : ! The determinants have now been reduced r — s orders lower 

■<3V 

ITllian at first, and the row of the denominator which is missing in 
l^Uhe numerator is that in which every constituent is unity. 

!|! Taking out the factors in the-numerator which occur through - 
I^Jut any column, we get 

loo 1 
H1 

/ I ....... I X 

r . ... . . . . r + s — I 

, . I 


= (-1)' 


,rt» r "" li .;/ + s-l)w-" j 

1 r I r — s ~ ’ ~ 7 

L — 1 -— ' 1....1 


r + 1 


r + s + i 


O+1) 00 .(r-rS+IjC® 


The row now missing in the numerator is the last, and in the 
numerator r occurs instead of r-f 1. 

Reduce the order of the determinants in the same way as 
before, s times, we get 

j r + s 


Substituting this value for r R s hi (ii), we get for the value 
of r A r 




.| 2r r 1 | r -f fi 
r 


(I r | sf 


s (n + r+ 1 )tr+*+n 


(A) 


From this formula we can without much labour determine 
the coefficient, in any numerical example, of the highest factorial 
in the expression for y. 

The values of the remaining coefficients are most easily ob¬ 
tained by an indirect method, as follows : — 

Let r B p be the value of r A r when aA 5 is substituted for y x1 
then 


A 


C r o(-l) f+S 


2r+ 1 \r + s 


(j 7 p) 2 \r-s + i) (r+s 




Now by aid of the formula 

Ux Vx - Vx~i I'jb- buss-2 A- 2 v x + &e. } 

we have 
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n (s) L-$ (?£ — I jw-i) -A--—_ T + , . . 


p+l 


(P+ 0 (P + 2 ) 


! 5 


+ (- i Y *|§g* (n -p) (s -A 

I S — p j P + P + I 


_ {pi -f l)(/ :> +i , + 1 ) + &c. 


Hence may be written in tlie form 


\2T +i I r + s 


__ l __ y P= s I- \ n - P l£ (n+ i)( p+p+V (_ i)r+» +* • 

(jr is ) 2 jr — 5 (» + '/*+ i)( r+ * +1 ) r-o j S _^ [ w _g jp+p+i 


or 

'** SSir (_ I'jP + r 

^ 30 'rp= 0 V ‘ ' 


\ 2r + i 


|r + s \ n~P .\P \ n + 


(l r) 2 i n + r + i '| s ) r—s \s — p\p + p+i \ n — p — p 


or, changing the order of the summations, 


rBp-^agC-i)^ 


How 


sr + s . 


12f-j- I jW—_£> |p j/fr-f I 


ir + s 


(jr) 2 j n + r -f I jp-t-p+i | n — p—p \ s j r — s j s — p 


If 


is = eoeff. of x s in expansion of j r (l—x) r ~ 1 , 


and 


i (i— $y-p 

( —i) r± * ■-— - ■ ■■ ■ is — coeff. of x r ~ s in expansion of- 


i r — s is — p 


I r—p 


when s has any value from s = p to s= r ; therefore 

\r+s i 


^:=A-o r+5 


is = coelf. of x r in expansion of 

! r 


I s i r “ s I S ~P 




that is, it is = 


Hence 


j r-P 


\r | r+p \r+p 

or 


\ r ~P \p j P f“P j P 


\2r+ i |p m+ i | n —p | r+p 


M P = (- iy '==■- ....—, 

p {[r] 2 + r + i [p+j?+ I j n — p—p \r—p j p 

which may be written 


rB P = 


: 2r + I j n + i | p / E V’ 


3 \n + r + 


(r+ i) c -p+ r -D, 

•f I V E / 
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ooi 

Jrhere 

E' applies only to (r + i)f-/»+ r - 1 ) and E' (r + 1 )(-*+»-*) = (r + 2y- p+r - 1 \ 

Sind 

<? | 

Si E applies only to nW and E # = (n + I )W. 

Hence 

j 2r + 1 | n+ I p { A — A!' 

r-Bp = 


(l r) 3 | ?7 + r + 1 \ E 
2f + I ; >1 + I | P p=r LT ( ” 1 / 


^jp=0 


(j r) 3 | n + r + 1 p_u [P | ^ — p 


) »W (r + iy-P+r-i), 

^_ r )(p) A /r ~^ (r+ iy-p+^-D, 


= ^_JL_ p(« ^ (r+ D<—^ 

(l r ) 3 1 ffl + r + I p |J> j r — P [P Z 


I 2r + 1 | n — r /[p \ 2 p _ 


__ , V ,=r ■■..-. .— __ \ \ 

{nT7+r\~p-r \\r_J \ n-r-p+p ir + p-p+l \p \r-p 


| n+i 


Now 


n + I 


; n - r— P 4- p r + p - p + I 


is = coeff. of a ,r+ ^“ iJ+1 in expansion of (1 •*-a?) n+1 > 


and 


therefore 


_=_is = coeff. of xp in expansion of (1 + x ) r ; 

j p \r — p 




| n+ l [r 


j ?7 —r — p+ p \ r + p-p + I i£ y — V 
is = coeff. of af +p ' K in the expansion of 1 4-*»* +r+1 » and is therefore 


! n + r -f I 


r + p + I 77 — p ! 


Hence 


A 


|n —r \2r+ 1 f _p\* 

— I 5 


|r + p + I 177 — p ip —r \_r, 
{pO°-op 

(r + />+i)k -*0 i p —r 


(iv) 


Now the original equations for finding ,A r) r A r _j, &e. may be 
written 
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. A, 23jt'!x» + r A,.,2;:; 




+ rAo s!=o xW-KZl . • — r A t+ i *« + «) = 0; 


so that when 


r A r , r A;-„j . . . r At4_i 

have been determined, r A t may be obtained in the manner that 
r A r was, but with 

y x - r A r x {r) ~- r A t+ i a? cm) 

written for y x , 

Hence by (A) 


rA*« 


4=0 


(-l)*+* 


\ 2 t+l \t + S 1 

([* | sf \ t — s [n + t + i)o+‘+u 


X 


— r A p . ,- r A t+1 ^+1)} 


= eh — <B r . } .A r — . r A r _i--. . . — i^p . r A p —. . . Jh +1 . r A f+ i 

— eh Av • rAjv, 

4 ‘ 4 *- t+1 (t + 2>+ 1 )^-^ |j>-* rAp ‘ 

Suppose £ =s r — i, then by (v) 

n — r 4-1 o . 

j- \r —1 r—lAj «—1 ' ~~ 1 r*h« 


(▼) 


Similarly, if in (v) t be put equal to r— 2, we get 


. . n—r+ 2 . f » —r+l _ . 1 

-h-2 = r^2Ar^2 ■ (f-I ) 2 ^r-lA r «i---—- f 1 r A r ^ 


2r 


/ 


(»-r+2p {r( 2 )p A 
(2r- I)W !2 " r 


•2 A r - 


n—r+2 


■ (r-i) 


- (r-1) 2 + 


(*- 


2 )(2)| r (2)}2 

Am 


2 (2r)(‘ 2) 


.In the same way if t be put equal to« — 3 in (v), it will be found 
that the resulting equation may be reduced to the form 

i‘Ar-3 “ f-aA r .3 —(r—2) a r * 3 A r .« 


4 ( W -r + 3 p i(r-l)ffl[ a 

: 2 (2r — 2)® r “ l1 r “ ! 


(^-r +3 p){rff)p 
13 (2/*>^ 



We see then that as far as ,.A r _ 3 the value of r A s is of the 
form 
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(B) 


Let as assume that all the coefficients ,.A r _ b r A.,._ 2 . . . r A s 
are of this form; we will prove that r A^ must also be of this 
form. 

For by (y) 

\ __ \ ^ P =r + a 

,*=»* {n-s + i)(P-s^){ p (r-s^Y2 <f=r , , 

= ^ + p)o>-«+i> ■/>-.*+I *=* 1 j ' (2 t)«-P) ,t~v 1 1 

(for, by hypothesis, B holds for all valaes of from s to r) ; therefore 

*p~r ^t=r ; (n~s+ i)(t-g-n ){^-«+ 1 )} 2 


>A «~ l S “ lA5_1 2 ^ =s ^ l ^ P+t ( 2 tyt-l 0 \ t-p (s + pfP-s+O IP-5+ I 

or, changing the order of the sammations, 




(rc--a + i )(wi)pa-*-n)>2 
{2t)0~P^ \ t—p (s+p) { P~ s+l) \p~ $ + 1 


At; 


At* 


Now 


7 - ■>• + 1 

J-p _p—s + 1 


is = eoeff. of x l ~v in expansion of (i + a?)*-** 1 , 


and 


(-I)PH 


Hence 


\t+P 


I t — s 1 p + S 


is = coefF. of xP +s in expansion of (l + x) 


w+s-l■ 


->=$-1 \ L J ^ _ 


t — S +1 


t+p 


t—p \P~ s + I t — s I p + 8 


is equal to the coefficient of x t+s in the expansion of (i+®)° 
and is therefore zero. Hence 


2;=a-ir 5 


or, finally. 


t+p 


t A — I 


I t — p\p — S+I\p + S \t-* + I 2.s — I 


,.A^ =,_ 1 A.: 1 -^ = ; ( -r)— < »-«+. .. mm , A , 


\21 


h~5-f I [28—1 


~ ( (aip^^ 7A-s~+ i“'~ cAs ’ 


which is of the form (B) with s—i in place of s. 
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■cTl 1 

loo i 
i^fi 

Bat r A r _ u r A r _ 2 have been already shown to be of the form 
;n{B), so therefore is r A r _ 3 , and so then ,.A r _* &c. Hence, for all 
lvalues of s from s = o to s = r — i, r A s is of the form given 

#y (B). 

The result of the foregoing investigation may be thus stated: 
! 3 f r A s be the coefficient of $ (s) obtained, by the method of least 
squares, from the n -{-1 values of y x when x = o, i, 2, . . . n, n 
being greater than r, and y x being of the form 



r A r 

+ x 0 ~~ l) 

. • • + -,A, A‘ ,s| + . . . + r A 0 ; 


then 






A. = 2 ( ‘= r ( 

j y+s ('*-«)“- 

' (2 t)«~ 

‘ A ‘’ 

(B) 

where 







, . j2 t+l 
- iV+* 

di ip) 2 ( n 

( ' + p)W 

+ ^Ip+ 1 > VX ' 

(A) 


The following are the values of the constants, calculated 
from these formulas, as far as they are required if fifth differences 
are assumed to be constant. 
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Constants in an Equation for y x . 


ijnupp. 1879. 


5 ° 


-> 

J 


loo 1 

HI If fourth differences are assumed constant, 5 A 5 must be put 
'^equal to zero; if third differences are taken as constant, 4 A 4 and 
;^;A 5 must be put equal to. zero ; and so on. 

;|; In employing these formula?, it will generally save labour if 
|lwe subtract from each of the quantities y the quantity given by 
;^he formula 


Cr .r 1 '' + C r _. u L .1 + C i,, 


where Ck, C ; ,_ l5 .... O 0 are so chosen as to make 


y x - (C, + C,-! ,r°'" 1 + ... - < id = o 

for all values of x from ce = n — r to x = n. Thus all the 
values of y x for which the coefficient is large will vanish, and 
the other values will also, as a rule, be very much diminished. 
The formulae (A) and (B) may then be employed, using the 
remainders, after subtraction, instead of the original values of 
y x ; and if r A’ s is one of the coefficients so determined, the 
coefficient in the original value of y x will be 

C, + j-Af 

The above investigation may also be modified, so as to 
apply to the case when the coefficients in the value for ?/,. are 
known to he connected by one or more linear relations, and 
also for approximating to the values when the relations are 
not linear. 

For, suppose it known that 

u r . ».A>. 4- o r —1 + . . . -f . ,.A p + . . . 4 - eta • rAo 4 - q = o, (vi) 

we must make 

2 Wo (»'A r SL'( r ) + r A,—1 ak + . . . + r A p xW + . . . + r A 0 —y x ) 2 

— A (Cty . r A -t- C r _ 1 . ,.A»— 1 + * * • + Up . r A p + . . • + Uq , r A 0 + q') 

a minimum. 

The resulting equations are of the same form as before, 
except that now 

A + 5 t . =0 ' t/ x 

occurs instead of 

^CC=() ^' P l!x * 

p having any value from o to r. 

Replacing 

KZl XP V* -hy A a p + 2*1? y x 
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in the solution, we get 


(■ 


■ i)p+* 


I (^ + p)(2p' 


([# |p ) 2 (» + t +1 )^+-°+ 1 ) 


{A^ + ^o &*y*}\ (vii) 


whence r A s &c. may be determined in terms of A. ^The value of 
A must then be determined, so that (vi) is satisfied ; and this 
value must then be substituted in the expressions for r A s , &c. 

If the relation be not linear, suppose <j) ( r A r , r A r _ 1? &c.) be 
the relation. 

Let r A' r , r A'._ 1? &c. be the values of the constants found by 
means of (A) and (B) without reference to the given relation. 
Suppose r A r = r A ' r -f b r A’ r &c., then 


(p (Ij.Aj.~J" 5 j-Aj.j yA-r—i 3" ® )-Aj-_xj &c.) — o 5 

therefore 

0 ( r Ag &C.) + 5 r M r -1" 5 j.A,._i -—-r &c. = o, 

d r A' r d } .A',._i 

or 


+ r A r -1 + . . . + 0 (.At, .At.!, &c.)-^At 44 -.At_! - &c. 

d r A‘ r d r Aj._i d r A‘ r cl r A‘ r ~.\ 


If we suppose that d,,A r &c. are so small that their squares 

and products may be neglected, then &c. maybe treated as 

constants ; the relation will now be of the form (vi) and we 
may find as above the values of r A r , r A r _ l5 &c. 

Meteorological Office, Toronto , Canada, 

1879 , April 30 . 


Royal Observatory , Greenwich, 
1879 ,'dtdy 29 . 

It is known generally to the Society that meridional observa¬ 
tions on a comprehensive scale are carried on at the Royal 
Observatory with unceasing regularity; but it is not so generally 
known that a considerable number of extrameridional or irre¬ 
gular observations have accumulated during my superintendence 
of the Observatory . 

I have thought, therefore, that I might with propriety lay 
before the Society an arranged list of observations of the last- 
mentioned classes, nearly the whole of which are included in the 
annual volumes of Greenwich Observations. 

£}. B. Airy. 


o, 
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